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Direct mass measurements of Cd isotopes show strong shell gap at N=82 
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Abstract 

A 238 U projectile beam was used to create cadmium isotopes via abrasion-fission at 410 MeV/u in a beryllium target at 
the entrance of the in-flight separator FRS at GSI. The fission fragments were separated with the FRS and injected into 
the isochronous storage ring ESR for mass measurements. The Isochronous Mass Spectrometry (IMS) was performed 
under two different experimental conditions, with and without Bp-tagging at the dispersive central focal plane of the 
FRS. In the experiment with Bp-tagging the magnetic rigidity of the injected fragments was determined by an accuracy 
of 2 • 10“ 4 . A new method of data analysis, using a correlation matrix for the combined data set from both experiments, 
has provided mass values for 25 different isotopes for the first time. The high selectivity and sensitivity of the experiment 
and analysis has given access even to rare isotopes detected with a few atoms per week. In this letter we present for the 
129,130,1310^ isotopes mass values directly measured for the first time. The Cd results clearly show a very pronounced 
shell effect at N=82 which is in agreement with the conclusion from y-ray spectroscopy of 130 Cd and confirms the 
assumptions of modern shell-model calculations. 
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1. Introduction 

Accurate mass measurements reflect details of the evolu¬ 
tion of nuclear structure and stability as well as the energy 
levels and spatial distributions of the bound nucleons [lj. 
A first microscopic explanation of the observed shell struc¬ 
ture and the corresponding magic numbers i! of neu¬ 
trons and protons, at which the nuclei have larger binding 
energies, refined the understanding of nuclear properties 
and manifested the application of the shell model. The ad¬ 
vent and application of radioactive nuclear beam facilities 
[d| and novel mass spectrometers Q have constantly en¬ 
larged the number of known isotopes with unusual proton- 
to-neutron ratios and thus explored the properties at the 
outskirts of the chart of nuclides and soon it became ev¬ 
ident that the nuclear shell structure can change towards 


* Corresponding author 

**Part of doctoral thesis at JLU GieBen (2015) 


the driplines. Shell quenching, complete disappearance, 
or even new magic numbers have been theoretically pre¬ 
dicted 0| and observed in experiments @, 0, 10, [Ill, [HI • 
The best known examples, both theoretically and experi- 
mentalhc are the N = 20 and N = 28 “islands of inver¬ 
sion” 0, El 3, El 0, E EH where the gain in correla¬ 
tion energy driven by quadrupole deformation is able to 
overcome the normal level population deduced from the 
standard spherical mean field. As a result the traditional 
N = 20 and N = 28 shell closures completely disappear. 
It has also been argued that such a shell quenching would 
occur for neutron-rich N=82 nuclei. This phenomenon 
was originally suggested based on Skyrme Hartree-Fock- 
Bogouliuvov calculations [6]. Almost simultaneously it 
was realized that the occurrence of an abundance trough 
around A ~ 115 in r-process calculations 0 could be 
cured using a mass model with a quenched shell gap far 
from stability 10 10, 20- The abundance trough around 
A ~ 115 is associated with a ‘saddle point behavior’ seen 
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in the two-neutron separation energies for Z « 40 and 
N = 75-82 in several mass models related to a transition 
from deformed nuclei around N ~ 75 to spherical nuclei 


at N=82 [21]. In mass models with a quenched shell-gap 
like the quenched extended Thomas-Fermi model (ETFSI- 
Q) [IS] the deformation is greatly reduced and conse¬ 
quently the ‘saddle point behavior’ in the two-neutron sep¬ 
aration energies disappears. However, it should be pointed 
that the ‘saddle point behavior’ and the quenching of the 
shell gap are not necessarily related 21 ] as the first could 
also be associated to instabilities of mean-field models in 
regions of shape coexistence that require the inclusion of 
beyond-mean field correlations [0. It is interesting to 
notice that the ‘saddle point behavior’ is absent in the 
Duflo-Zuker shell-model inspired mass model [23] and con¬ 
sequent does not result in an abundance trough around 
A ~ 115 despite of predicting a relatively large N=82 shell 
gap far from stability. Furthermore, the strong abundance 
troughs predicted in calculations were largely due to the 
use of the so called classical r-process approximation which 
assumes an ( n , 7 ) ( 7 , n) equilibrium and an instanta¬ 
neous freeze-out of neutron captures E3- Nowadays, we 
know that the r-process occurs in a very dynamical envi¬ 
ronment with continuously changing astrophysical condi¬ 
tions and that the competition between neutron captures, 
beta-decays and fission has an important impact in the 
final abundances [230. 


Given the relevance to both r-process nucleosynthesis 
and nuclear structure of the existence of a new ‘island of 
inversion’ approaching Z ~ 40 for N=82, there have been 
many experimental attempts to provide evidence for the 
quenching of the shell gap. It should be clarified that 
quenching is sometimes used as a synonymous of reduction 
of the shell gap. However, in the context of this letter 
we will use ‘quenching’ according to the most common 
usage of disappearance of shell effects pointing to an ‘island 
of inversion’ analogous to the well-studied cases around 
32 Mg for N = 20 and 42 Si for N = 28. The experimental 
information about the evolution of the N=82 gap is so far 
mainly indirect. 

The present Isochronous Mass Spectrometry (IMS) ex¬ 
periment and the previous Penning trap mass measure¬ 
ments for the tin isotopes f 2 fl yield direct information on 
the shell effect. We have applied IMS to determine the 
masses of 129 d30,i3i(^ [ j i so topes. The cadmium and tin 
results provide a direct experimental evidence for the evo¬ 
lution of the N=82 shell gap south of 132 Sn. 

2. Experiment and data analysis 

A 410 MeV/u 238 U projectile beam was extracted from 
the synchrotron SIS-18 27] with an average intensity of 
lT0 9 /spill and impinged on a 1032 mg/cm 2 beryllium tar¬ 
get at the entrance of the fragment separator FRS (28| . 
Neutron-rich fission fragments created via abrasion-fission 
were separated in flight with the FRS applying pure mag¬ 
netic rigidity (Bp) separation with the standard ion-optical 


operation mode. The separation mode, without degraders, 
was enabled by the large mean velocity difference of the 
projectile fragments and fission products and the restricted 
angular acceptance of the FRS. Practically this means, a 
suitable -Bp-selection with the FRS can provide fission- 
fragment beams without significant contributions of pro¬ 
jectile fragments. The ions of interest were injected into 
the Experimental Storage Ring ESR [0 for IMS |3p], 0 
at a mean velocity corresponding to the transition energy 
of 7 1 — 1.41. The magnetic fields of the FRS and ESR 
were set for 13 3,i35,i36g n i ons j n diff eren t runs, i.e., these 
isotopes were subsequently centered at the optical axis. 

The ESR was operated in the isochronous mode [0 
without application of any cooling. This means that the 
velocity spread of the fragments is determined by the Bp 
acceptance of the ion-optical system. In a previous publi¬ 
cation 33] we have demonstrated that for IMS experiments 
in addition to the revolution time of the stored ions the 
magnetic rigidity or velocity measurement is required, be¬ 
cause the isochronicity is strictly realized only for a single 
mass-over-charge ( m/q ) value. 

In principle, this additional measurement is not required 
for Schottky Mass Spectrometry (SMS) because the rela¬ 
tive velocity spread of the different stored and cooled ions 
can be as low as 10~ 7 . Nevertheless, our refined SMS anal¬ 
ysis has revealed that an additional influence of the mean 
velocity causes an observed correlation Q. 

The method of IMS including Bp-tagging can be illus¬ 
trated by the simple first-order formula 


d(m/q) _ 2 d(T) 

-7— “ 7 —HI- 2 

m/q T 7 1 
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where T, 7 t , and 7 are the revolution time, the transition 
energy, and the relativistic Lorentz factor, respectively. 

The additional velocity (v) and magnetic-rigidity mea¬ 
surements in FRS-ESR IMS experiments require special 
methods due to the ESR operation with fast extracted ion 
bunches characterized by a width of (0.2 — 0.5) ps. Par¬ 
ticle detectors inside the FRS would have severe problems 
to identify event-by-event the fragments and accurately 
measure v and Bp. A Bp-resolution of 10 ~ 4 or better is 
required to achieve a mass resolution of about 200 keV for 
m/q close to ideal isochronicity, details are presented in 


reference [35|. In this context, one has to take into ac¬ 


count that the FRS transmission is d{Bp)/(Bp) = 2% and 
the corresponding ESR injection acceptance is roughly one 
order of magnitude less. Therefore, mechanical slits with 
a slit opening of ±0.5 mm placed at the central dispersive 
focal plane of the FRS were used in a pilot IMS exper¬ 
iment 0]. The slits defined in this way the magnetic 
rigidity (Bp-tagging) of each injected ion with an accu¬ 
racy of 2 • 10~ 4 . The operation of the new isochronous 
Rare RI-Ring at RIKEN 0] can easily implement addi¬ 
tional Bp and v measurements event-by-event, because of 
the DC beam from the cyclotron. IMS measurements have 
been performed with and without Bp-tagging for the same 
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settings of the magnetic fields of the FRS and ESR. The 
revolution time of the circulating ions in the ring were 
measured with a time-of-flight (ToF) detector equipped 
with a thin carbon foil and two micro-channel-plate (MCP) 
branches [32J placed in a homogeneous magnetic dipole 
field of about 8.4 mT. The created secondary electrons 
in the carbon foil were isochronously deflected onto the 
MCPs to generate timing signals at each turn. The signals 
were recorded with commercial digital oscilloscopes (Tek¬ 
tronix TDS 6154C, 40 GS/s, 15 GHz; LeCroy LC584AM, 
4 GS/s, 1 GHz) 

The data sets of the two different experiments, with and 
without Bp-tagging, were combined and analysed with 
the new correlation-matrix method [Hil. In this work 
we used for the first time a variable correlation factor 
dependent on the measured y . A detailed description of 
this novel analysis will be presented in a forthcoming pub¬ 
lication g0i. The advantage of this new analysis is that we 
could include ions with very low statistics down to a few 
events of a single isotope. In this way we can now present 
more than 20 new mass values which were not included in 
previous IMS evaluations of the same experiments 4l|, [36]. 


3. Results and discussion 


The first check of the reliability of the new data anal¬ 
ysis is the determination of the systematic error of the 
combined experiments. An average of 1600 events per ref¬ 
erence mass could be applied in the analysis with the com¬ 
bined data sets. We determined the systematic error by 
using all reference masses and subsequently treat each of 
these masses as being unknown. The procedure can be 
illustrated by equation [ 2 ] and figure Q] 


n 
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where are the mass values and a( e ^ the uncertainties 

of the reference nuclides. af tat are the statistical errors of 
the measured masses Mj. N n is the number of reference 
masses and cr syst the systematic error. In this analysis 47 
reference masses have been used Phi . The investigation 
covers the full m/q range of the present analysis. The 
mean value of the projected distribution is 1.29 keV. The 
deduced systematic error is 172 keV (standard deviation). 
For most of the new mass values the systematic error 
is the dominant contribution to the total error which 
results from the sum of the variances. The statistical 
error is determined by the uncertainty of the revolution 
time, which strongly depends on the number of turns 
for which the timestamps have been recorded. In the 
matrix method the statistical error for each nuclide was 
calculated by the square root of the diagonal elements of 
the inverse matrix. Details will be presented in reference 
0 . 



m/q 

Figure 1: The systematic error of our mass measurement was deter¬ 
mined by reanalysis of each reference mass, i.e., each reference mass 
is treated subsequently as an unknown species. The distribution of 
this analysis is depicted as a function of mass-over-charge values and 
shown in the insert as a projected histogram. 


A main motiation of the present mass measurements is 
the investigation the N=82 shell gap for Cd isotopes. 
We briefly review previous experiments and conclusions 
before presenting our results. Ref. 43] determined the 
beta-decay half-live of the 130 Cd nuclide and its Qp value. 
The mass of 130 Cd in the Atomic Mass Evaluation 2012 
(AME12) gj is based on this Qp value. It was concluded 
that the large measured Q^-value of 8.32 MeV is an 
indication of shell quenching because the experimental 
value was in good agreement with the predictions of the 
ETFSI-Q mass model gjj. However, it is about 1 MeV 
larger than predicted by other mass models such as the 
Finite Range Droplet Model gl or the Duflo-Zuker 
model 231 ] . A recent experiment [47] has determined 
a substantially shorter beta-decay half-live with more 
statistics for the 130 Cd nuclides, which questions the 
determination of the Q^-value in [431 ]. The measurement 
of the 2 + excitation energy for 128 Cd showed that their 
value was the smallest of the known N=80 isotones [481. 
Again, it was interpreted as evidence of shell quenching. 
However, beyond-mean field calculations g)j showed that 
such a low excitation energy could be explained by the 
presence of prolate configurations in the wave function of 
128 Cd. This may explain the measured small magnetic 
moment observed in 127 Cd [50] . The same calculations 
predicted a much larger excitation energy of the 2 + level 
in 130 Cd. This excitation energy has been measured 
to be 1.321 MeV [5l| in agreement with shell-model 
calculations that assumed a N=82 shell closure. This fact 
together with the large excitation energy of the 2 + state in 
126 Pd [52j (1.311 MeV) has been interpreted as evidence 
for a robust N=82 shell closure for both cadmiun and pal¬ 
ladium isotopes. Recently, a 7 -spectroscopy experiment 
53} has determined the 3/2“ state in 132 In corresponding 
to the p 3 / 2 proton-hole state in 132 Sn. It is interesting to 
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Table 1: Fitted liquid-drop parameters. The parameters, in units of 
MeV, for the Weizsacker formula have been obtained by a fit to the 
mass values of reference t44l . 



Figure 2: Experimental shell effect. Difference of the measured mass 
values and the smooth Weizsacker formula given by Eq. [3] for the 
elements Sn and Cd. The parameters used in the formula are pre¬ 
sented in Tabled The curves clearly show the extra binding energy 
due to the contribution of shell structure for the Sn and Cd nuclei 
in this mass region of N=82. 


note that these new data are in perfect agreement with the 
predictions of the Duflo-Zuker monopolc hamiltonian [55j 
which is the base for the Duflo-Zuker mass model. The 
new experimental data has allowed to determine all 
relevant proton single particle states showing that no 
proton subshell closures are expected for light N=82 
isotones pointing to a significant reduction of the N=82 
shell gap in analogy to the reduction observed in isotones 
heavier than N=82. This reduction does not necessarily 
mean a quenching of the shell gap. It is well known that 
shell gaps reach a maximum value for magic numbers of 
protons and neutrons [561. 


Figure [5] shows the difference of the experimental mass 
values for Sn and Cd isotopes and the prediction of the 
liquid drop model 1, [57|]- The liquid drop parameters 
have been deduced from a fit to the tabulated values of 
the Atomic Mass Evaluation 2012 (AME12) [i3 |. 

The liquid-drop binding energy Bld has the form 


Bld 


b vo iA - b surf A 2 / 3 


1 L {N~Z) 2 

2 bsym A 


3 Z 2 e 2 
5 R c 


0 , 

bpairA _1/2 , 

O . U . A- 1/2 
" u pair- ry i 


e-e nuclei 
e-o or o-e nuclei 
o-o nuclei 


(3) 


with R c = 1.24 fm -A 1 / 3 . 


In figure [2] the differences between experimental mass 
excess (ME) and the corresponding liquid-drop values are 
plotted. The ME is the difference between the atomic 
mass and the corresponding mass number, both expressed 
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in atomic mass units. The liquid-drop values are calcu¬ 
lated with equation [3] and the parameters in table □ The 
strong extra binding energy for the neutron-rich Cd iso¬ 
topes is clearly observed at N=82, it amounts to roughly 
-14 MeV. The error bars for 129 430 qj are we n w ithl n the 
symbol size whereas for 131 Cd we have a large statisti¬ 
cal error because only two events have been recorded with 
few number of turns. The experimental results in figure 
[2] clearly demonstrate that the tin and cadmium isotopes 
are characterized by the same difference of shell correc¬ 
tions at N=82. This observation contradicts reference jidj 
and clarifies the situation of the different statements in 
references [HI, [HM E3]. 

The experimental ME values for 129 . 1 30,131 q c j nuc i e j anc j 
the errors are presented in table [2] Our directly measured 
mass value for 130 Cd atoms is 600 keV lower than the value 
deduced from Qp measurements [43}. This difference could 
be caused by a contribution of an excited state in the (3- 
decay experiment. In principle, also in our measurements 
for 129,131 Cd nuclei unresolved isomers could contribute 
fsdl . [Hd} whereas in 130 Cd the lifetime is too short to inter¬ 
fere in this experiment. However, in this case the observed 
strong shell effect for Cd isotopes at N=82 would be even 
stronger for pure ground-state masses. The nearly equal 
shell effect for Sn and Cd is a striking feature of the present 
experimental results. 

Other established results from mass measurements to 
discuss the evolution of nuclear structure and shells are the 
nucleon separation energies and their derivatives. The one- 
neutron separation energy S n and the pairing gap energy 
is defined by: 


S n = -M(Z, N ) + M(Z, N - 1) + Af (0,1) (4) 

Ai 3) = (-1)^ [S n (Z, N ) - S n (Z, N + 1))] (5) 


A and Z are the mass and proton numbers, respectively. 

In the upper panel of figure [3] the shell evolution is man¬ 
ifested by our experimental (indicated in red) or AME12 
one-neutron separation energies and their derivative near 
N=82. In the lower panel, we display an estimation of the 
gap representing the single particle spectrum (3-mass for¬ 
mula) which in principle is the derivative of the S n -values. 
Both presentations demonstrate a strong shell effect at 
N=82 for Sn and Cd isotopes in complete agreement with 
the characteristics of figure [2] Again it is a striking exper¬ 
imental observation that Sn and Cd isotopes are governed 
by the same strong shell gap at N=82. The magnitude of 
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Figure 3: In the upper panel experimental one-neutron separation 
energies ( S n ) are shown for Cd and Sn isotopes. In the lower panel 
the corresponding pairing-gap energies are presented. 


Table 2: Measured mass excess values (ME) of Cd isotopes. The 
statistical (cr s tat) and the total (cr tota i) errors are tabulated. The 
systematic error is 172 keV for all the new masses measured. In ad¬ 
dition, the extrapolated (jf) values [3 and the 130 Cd mass deduced 
from Q^-measurement [4311 are listed. 


isotope 

ME 

[keV] 

& stat 

[keV] 

&total 

[keV] 

ME AM Eli 

[keV] 

counts 

i2y Cd 

-63145 

17 

173 

-63509(196)# 

18 

130 Cd 

-62131 

41 

177 

-61534(164) 

5 

131 Cd 

-55583 

925 

941 

-55331(196)# 

2 


the shell gap is within the experimental uncertainties the 
same. 


A comparison of the experimental values for Cd isotopes 
with different theoretical mass models is presented in fig¬ 
ure [U The models are based on microscopic-macroscopic 
descriptions HE EE EE EE , the Hartree-Fock-Bogoliubov 
(HFB) theory ||6lj, and the shell-model inspired model of 
Duflo-Zuker 



Figure 4: Comparison of experimental mass values for Cd with dif¬ 
ferent theoretical predictions 


In general, the comparison with theoretical predictions 
clearly demonstrates that the agreement with experimen¬ 
tal results is getting worse towards the shell closure at 
N=82 and of course in the n-rich region, where no experi¬ 
mental results were available before. There are basically 3 
different regions in this comparison. In the first region for 
N<81 the theories are relatively similar in their trend but 
all predict too low mass values (over-binding). At the shell 
the predictions are quite unreliable and beyond N> (82), 
where no experimental data exists, they widely scatter. 
The upper panel show the described comparison for Cd 
isotopes and the lower panel the corresponding results for 
Sn isotopes. The masses of the presented Sn isotopes have 
been measured with Penning traps with higher resolving 
power and should not have influenced by unresolved iso¬ 
meric states. The predictive power of the different mod¬ 
els can be quantitatively characterized by the a rms values. 
They are listed in table [3] for Sn and Cd isotopes in the ex¬ 
perimentally known range of AME12 and for Cd isotopes 
presented in figure [4] 

The comparison with theories over a larger Z range for 
N = 82 isotones is presented in figure[5]for the one-neutron 
gap, Ai n (Z, 82) = S n (Z, 82 ) — S n (Z, 83). Experimental S n 
values deduced from mass measurements exist only above 
cadmium. The figure clearly shows the large gap measured 
for 130 Cd that is within the experimental uncertainties of 
similar magnitude than the one for 132 Sn. No theoretical 
model is able to reproduce such a large gap for Cd and Sn 
isotopes. Looking at the evolution of the one-neutron gap 
and averaging the odd-even staggering one observes that 
only the ETFSI-Q and HFB27 models predict a strong 
reduction towards Z=40. 


The new data for Cadmium isotopes is not consistent 
with an early quenching of the N = 82 shell gap as claimed 
before by several publications iE 481. However, it cannot 
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Table 3: Comparison of measured data with models. RMS devi¬ 
ations for Sn and Cd isotopes for different theoretical models are 
presented. In this comparison the new IMS values and the tabulated 
experimental values of AME12 are included. 0y mSj i,<7 rmSj 2 represent 
the values for the cases of all experimentally known Cd isotopes and 
the Cd isotopes in the mass range as shown in figure 4, respectively. 



[keV\ 

&rms, 2 

[keV] 

FRDM [46] 

619 

1537 

HFB-27 [61] 

523 

775 

DZ28 [23] 

315 

958 

ETFSI-Q [45] 

342 

512 

ETFSI-1 [60] 

479 

1200 

WS3 [59] 

396 

709 


nounced shell effect at N=82 which is assumed by modern 
shell-model calculations. The experimental values for the 
shell corrections are roughly the same for Sn and Cd iso¬ 
topes. 

The goal of future IMS experiments will be to measure 
the lower elements (40 < Z < 48) at and near the N=82 
shell. Only the mass values of the lower elements will 
be decisive for the observation of possible shell quenching 
as it is proposed in several theoretical models. However, 
experimentally this interesting region can only be accessed 
by the next generation exotic nuclear beam facilities [62 . 

IMil El. 
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Figure 5: One-neutron gap (see text) for N = 82 isotones. 


rule out a disappearance of the shell gap for lighter iso¬ 
tones. Therefore, it will be necessary to have experimental 
access to isotones with Z < 42 for which the differences 
between models become substantial. 

4. Summary and Outlook 

IMS experiments have been performed with and with¬ 
out Bp-tagging at the FRS-ESR facilities at GSI. A new 
method of data analysis using the correlation matrix for 
the combined data of both types of experiments provided 
25 new mass values in range of Ge to Ce (A=86-154) even 
for isotopes measured with very low statistics of a few 
atoms per week. The latter condition clearly demonstrates 
the high sensitivity and selectivity of the experimental 
method. In this letter we have presented the masses of 
the 129 > 1 30,131 q ( j i so topes which have been directly mea¬ 
sured for the first time. The goal was to investigate the 
evolution of the shell gap of Cd isotopes compared to Sn 
isotopes at N=82. The Cd results clearly show a very pro- 
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